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ABSTRACT 

By coupling a solar surface flux transport model with an extrapolation of the he- 
liospheric field, we simulate the evolution of the Sun's open magnetic flux and the 
heliospheric current sheet (HCS) based on observational data of sunspot groups since 
1976. The results are consistent with measurements of the interplanetary magnetic field 
near Earth and with the tilt angle of the HCS as derived from extrapolation of the ob- 
served solar surface field. This opens the possibility for an improved reconstruction of 
the Sun's open flux and the HCS into the past on the basis of empirical sunspot data. 

Subject headings: solar-terrestrial relations - Sun: activity - Sun: magnetic fields 
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1. Introduction 

The Sun's open magnetic flux is the part of its flux which is not contained in closed loops, but 
extends into the heliosphere. It is the source of the hel iospheric rnagne tic field (HMF) whose varia- 
tions are an important source of geomagnetic act ivity (lPulkkinenll2007l l and control the production 
of cosmogenic isotopes by galactic cosmic rays (iBeerlbonnl ^TA crucial feature of the HMF is the 
heliospheric current sheet (HCS), the interface separating the opposite polarities of the HMF. The 
tilt angle of the HCS (defined as the mean of the maximum northern and southern extensions of 
the HCS) is a key parameter for the modulation o f the flu x of galactic cosmic ra.ys in the inner he- 



liosp here (jKota Sz Jokipiilll983l : lFerreira fc Potgieter,i2004 lAlanko-Huotari et al.ll2007l : iHeber et al 



At a give n distance from the Sun, the HMF has an almost uniform magnitude in latitude 
and longitude (jBalogh et al.lll995l ). Therefore, its radial component near Earth, w hich has been 



measured by spacecraft sinc e the 1960s, faithfully represents the Sun's total open flux (| Owens et al 



20081 : iLockwood et al.ll2009l ). The tilt angle of the HCS could, in principle, be measured by multiple 
spacecraft orbiting at different heliolatitudes. However, with the exception of the Ulysses probe, all 
measurements of the HMF have been obtained near the ecliptic plane, so that direct measurements 
of the tilt angle of the HCS are not available most of the time. Therefore, such data are derived 
by extrapolation of solar surface field maps, such as those taken at the Wilcox Solar Observatory 
(WSO) since 1976. This yields the current sheet distribution at the source surface, where the field 
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is assumed to become radial (e.g.. iHoeksema et al.lll982l ). In the inner heliosphere, we may ignore 
dynamical effects such as the acceleration of slow plasma and the d eceleration of fast plasma that 
occur when neighboring parcels of plasma interact (jRilev et al.ll2002l ). so the field lines are assumed 
to stay purely radial beyond the source surface. Under these conditions, the morphology of the 
HCS may be inferred from the position of the current sheet at the source surface. 

A limitation of this semi-empirical determination of the HCS tilt angle arises from the de- 
creasing reliability of the surface field measurements at higher latitudes and from magnetographic 
saturation effects, which cannot be corrected without further assumptions. Flux transport mod- 
els based upon observed large-scale magnetic flux emergence (e.g., in sunspot groups) provide a 
complementary possibility to obtain information about the high-latitude surfa ce fields, which con- 



trol the open fiux and th e HCS tilt angle during solar minimum periods (e.g. IWang et al. 



1989a 



2000 



Mackay et al.ll2002l ). ISchiissler &: BaumannI (j2006l ) showed that such an approach reproduces 



well the HMF over multiple solar cycles, provided that the heliospheric curren t sheet is explicitly 
included in the field extrapolation (current sheet source surface method, cf. IZhao &: Hoeksema 
19953) 



The present paper serves a two-fold purpose. Firstly, it extends the study of lSchiissler &: Baumann 



(1200a) until the current solar minimum period and includes an explicit account of the HCS and its 
tilt angle, which can be compared with the corresponding observational quantities. Secondly, such 
comparison provides the validation for the application of these methods to reconstruct the open 
flux and HCS tilt angle backward in time. 

The paper is organized as follows. In Sect. 2, we give the description of the models used. The 
results are presented in Sect. 3: photospheric flux distributions in Sect. 3.1, the solar open flux in 
Sect. 3.2, and the heliospheric current sheet in Sect. 3.3. We give our conclusions in Sect. 4. 



Methods 



2.1. Surface flux transport model 



The surface flux transport (SFT) model describes the evolution of the magnetic flux distribu- 
tion at the solar surface as a combined result of the emergence of bipolar magnetic regions (BMRs), 
flux c ancellation, and transport by surface flows. The evolutioii of th e radially orientated surface 
fleld (|Wang k Sheelevlll992l : ISolankilEggj : IPetrie k Patrikeeva Ibood ) is controlled by latitudinal 



differential rotation an d meridional flow, t ogether with turbul ent diffusion due to granulation and 
supergranulation (e.g.. lWang et al.lll989bl : iMackav et al.ll2002l ). 



The SFT uses the radial component of the induction equation in the form 
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where S{9, (p, t) is the source term describing the emergence of new mag netic flux and Dr('nrB ) 
is the dec ay term para i neter izing the radial difl'usion of the magnetic field (jBaumann et alJ l2006l ) . 
Following iJiang et al.1 (|2009l ). the horizontal diffusivity r//j and radial diffusivity rjr are taken as 
600 km^ s~^ and 100 km^ s~^, respectively. For the meridional flow v{9), we adopt the profile 



-vo sin(2.4 * (90° - 9)) 15° < 9 < 165° 
otherwise, 



(2) 



wher e = 11 m s ^. This profile is largely consistent with helioseismic results ( Gizon &: Duvall 
200411. Fo r the l atitudinal differential rotation uj{9), we use the empirical profile determined by 



Snodgrasd (|l983l ). 



The basis for magnetic flux input to the SFT model by emerging BMRs is the USAF/NOAA 
sunspot group record Q|. Since a sunspot group typically appears more than once in the record, we 
consider a group only at the day of its maximum area. Figure 1 gives the monthly number and the 
latitude distribution (butterfly diagram) of the emerging sunspot groups for the time interval 1976 
- 2009, which provide the flux input for the SFT model. 

The observed sunspot areas have been conve rted to the areas of the corresponding BMRs 
following the procedure of iBaumann et al.l (|2004l ) . We include the magnetic flux contained in 
faculae and plages by er nploying the stat i stical relationship between sunspot area, Ag, and facular 
area, Aj, determined bv lChapman et al.l (jl997l ). 



Af = 414 + 21^ - 0.0036yl^ 



(3) 



(in units of millionths of the solar hemisphere) and take ^bmr = Ag + Af as area of the correspond- 
ing BMR. Since there is no information about the magnetic polarity in the USAF/NOAA data, 
we use Hale's polarity law to infer the polarities of the leading and following parts of the BMRs. 
We resolve the ambiguity arising from the overlap of cycles around activity minima by assuming 
that BMRs emerging below ±15° latitude during the overlap period belong to the old cycle while 
all others belong to the new cycle. The angular separation, A/? (in degrees), between the leading 
and following polarity patches of a BMR is assumed to be proportional to the square root of the 
BMR area (given in square degrees): A/3 = 0.6-v/Abmr. The angular separation is separated into 
latitudinal and longitudinal components, depending on the BMR tilt angle, 7, with respect to the 
E-W direction. We assume the relation 7 = 0.15A , which is consistent with observational results 



^http://solarscience. msfc.nasa.gov/greenwch.shtml 
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( Howardll99ll : ISivaraman et alJll999l . Dasi Espuig et al., in preparationjj. Finally, we calibrate the 
conversion factor between BMR area and magnetic flux by matching the observed and simulated 
values of the disk-averaged unsigned flux density, Bg = J J \ Br \ d0cos(A)dA/47r. 



2.2. Field extrapolation model 



In order to determine the coronal and heliospheric magnetic field from its source in the photo- 
spheric field, a field extrapolation method is required. For the field distribution on a global scale. 



the most widely used approach is the potential field source surface (PFSS) model (jSchatten et al 



196S 



Altschuler Newkirklll969l ). However, the PFSS model (which includes only volume currents 



beyond the spherical source surface, where the field is assumed to beco me purely radial) do es neither 
reproduce the latitude-independent radial field found with Ul ysses (IBalogh et al J 1 19951) no r does 
it match the measured interplanetary radial field near Earth (Schiissler &: Baumannl 120061 ) . The 



current sheet source surface (CSSS) model (|Zhao &: Hoeksema 



1995a 



Zhao et al 



2002), which 



explicitly takes into account the existence of the HCS, d oes not suffer from these de ficiencies and 
provides a reasonable match to the measured quantities (jSchiissler &: Baumannll2006l ). 



We briefly describe the main features of the CSSS model as follows. To include the effects of 
volume and sheet currents, the exterior of the Sun is divided into three parts, which are separated by 
two spherical surfaces, the cusp surface at r = Res, and the source surface at r = R^s (Res < R-ss)- 
In the region inside the cusp surface, the field is potential. In the region between Res and Rsst 
all flux loops are reconfigured with volume currents and current sheets, so that the field becomes 
completely open. In the region beyond Rgs, the field is purely radial. Apart from Res and Rss, 
the third adjustable pa,rame ter of the CSSS model is the height scale, a, of the horizontal current. 
Following IZhao et al.l (|2002l). we take a = 0.2R(r) in our calculations. We choose Rss = W-ORq 
according to the estimate of iMarsch &: Richter I (119841 ). The value of i?cs = 1-8Rq is determined by 



the comparison between the magnetic flux at the cusp surface and the measured radial field from 
0MNI2 data in Sect. 3.2. 

Technically, the opening of the flux beyond -Res and the introduction of the current sheet (s) 
is carried out by first calculating a global potential-field extrapolation for the entire volume above 
the solar surface. The unsigned radial component of the field at the cusp surface is then used as 
the bo undary condition for c alculating a magnetic field distr i bution between i^es and i?ss of the 
form in iBogdan &: Low! (119861 ) and used by IZhao &: Hoeksemal (|l995al ) . The orientation of the field 



Schiissler fc Baumann ( 20061 ) attempted to fix the latitude-dependence of the tilt angle independently by cali- 



brating with the longitude-averaged surface field integrated over latitude, viz. J \ J Brdcj> \ cos(A)dA/47r. However, 
some exploratory experiments we performed have shown that this quantity is rather sensitive to the observed scatter 
of the BMR tilt angles about the latitude-dependent mean. Since the effect of this scatter is not significant for the 
quantities that we study in this paper (open fiux and tilt of the heliospheric current sheet), we do not consider BMR 
tilt angle fiuctuations here. 
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lines is then changed where necessary so that the sign and magnitude of the radial field at Rc 
continuous. The distribution of the field outside Res can then be determined. 



IS 



3. Results 



3.1. Photospheric magnetic field distributions 

The initial condition for the photospheric flux distribution is the same as the one assumed 
by iBaumann et al.l (|2004l ). It satisfies an approximate balance bet ween the effects of poleward 
meridional flow and equatorial diffusion (Ivan Ballegooijen et al.lll99a ). The memo ry of the system 



regar ding the initial field depends on the value of the radial diffusion parameter, r]r (jBaumann et al 
2000 ). We use r/^ = 100 km^ s~^, which leads to a memory of about 20 yrs. We start all our 



simulations from the beginning of the USAF/NOAA sunspot group record in 1874, but consider 
the results only for the time period 1976-2009. Therefore, there is no remaining influence of the 
initial fields on the results presented below. 

The upper panel of Fig. 2 shows a comparison of observed and simulated time evolution of the 
averaged unsigned flux density at the solar surface, Bg (defined in Sec. 2.1), which has been used to 
calibrate the relation between area and magnetic flux of the BMRs providing the input for the flux 
transport model. Note that, without any other parameter adjustment, the model reproduces well 
the ratio between the maximum and minimum values as well as the very low surface flux during 
the current minimum. 

The lower panel of Fig. 2 gives the corresponding time evolution of the high-latitude surface 
field, averaged over the caps poleward of ±75° latitude. The field amplitude an d reversal times 
before 2002 are consistent with the observational results given by lArge et al.l (120021 ). For the current 
minimum, the reported values of the polar field do not give a consistent pictur e, reflecting the 
difficulties and limitations of polar field measurement. iPetrie Patrikeeva I (|2009l ) obtained polar 
fields of 5 - 6 G by analyzing photospheric and chromospheric vector polarimetric data obtained 
with SQLIS at NSO, w hich is consistent with our model. On the other h and, MWO magnetograph 
(jSvalgaard et al.l l2005l ) data and more indirect indices (jSchattenl |2005| ) suggest that the present 
polar field could possibly be a factor 2 smaller than that during the previous activity minimum 
around 1997 (see further discussion in Sect. 3.4). 



3.2. Open flux and near-Earth radial field 

We can calculate the total open flux ^open resulting from our extrapolation model by integrat- 
ing the unsigned radial magnetic field over the source surface, viz. 



^open(t) = Ri j j \BriRss, K <P)\dn. 



(4) 
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Actually, in the CSSS model the open flux is already fixed at the cusp surface, Res, smaller values 
of which leading to a larger amount of open flux. In order to compare with data provided by 
near-Earth measurements, we also calculate the longitudinally averaged unsigned radial fleld near 
Earth (at te = 1 AU and A = 0, thus ignoring the slight variation due to the angle of about 7 
degrees between the equatorial plane of the Sun and the ecliptic plane). 



BEit) 
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Rss 
TE 



2n 



\Br iRss,0,(f>,t) 



(5) 



Since extr apolations with the CSSS ni odel reproduce the latitude-independence of the radial field 
at r = te ( Schiissler &: Baumann 20od ) . we have ~ ^open/(47rr|) to a high degree of accuracy, 
so that Be also represents the open flux. Similarly, it has been shown that, possibly after some 
correction for kinematical effects due to the solar wind, the measured radial field near Earth can 



be taken as a reliable proxy for the total open fiux ([Owens et al.l 120081 : iLockwood et al.ll2009l ). 



Figure 3 shows 27-day averages of Be from our combined SFT/CSSS model (with Res = I.SRq, 
Rss = IORq, a = 0.2Rp), red curve) in comparison with the measured radial field from 0MNI2 data 
(blue curve) iLockwood et al. ( 2009 ) have suggested that the observed data better represent the 
open flux of the Sun if a correction of kinematical effects due to the longitudinal structure of the 
solar wind is applied. Data modified in this way are also shown in Fig. 3 (green curve). 

The phase relation between the solar activity cycle and the near-Earth field (and thus the open 
flux) is well re presented by our rn o del, with B^^ reaching its peak values ~ 2 - 3 yr after activity 



Mackav et al.l l2002l : IWang et al 



2OO2I : ISchiissler &: BaumannI 120061 ). Including the 



maximum (cf. 

diffusion in radial direction in the SET model and the realistic tilt angles of BMRs with respect to 
the E-W direction are the main reasons that lead to the correct phase relation. The only periods 
showing a significant disagreement between our model and the data are the ascending phases of the 
activity cycles, where the model values are too low. Possibly the model misses open flux from small 
coronal holes at intermediate latitude during this phase. In principle, this could be corrected by 
putting the cusp surface nearer to the Sun during this cycle phase or by assuming a non-spherical 
cusp surface. The amplitude of the variation of Be is also affected by the value of radial diffusion 
in the SET model. However, such parameter tuning would not provide further physical insight 
while the overall agreement between model and data is already encouraging and sufficient for the 
purposes of this paper. 



3.3. Heliospheric current sheet (HCS) 

Figure 4 shows the distributions of the magnetic field on the solar surface, on the cusp surface, 
and on the source surface during typical solar minimum and solar maximum conditions in our 



■^http: / /omniweb. gsfc.nasa.gov/ 
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coupled SFT/CSSS model. Around solar minimum (left panels) , the field at (and outward of) the 



source surface assumes a 'split-monopole' structure (IBanaszkiewicz et al.l 19981 ). the HCS separating 
the two polarities being located near the equatorial plane ( Hu et al. 20081 ). Near activity maximum 
(right panels), the HCS shows strong excursions in latitude and additional localized current sheets 
may occur. Apart from the current sheets, the field is always largely latitude-independent. The tilt 
angle of the HCS is defined by convention as the arithmetic mean of the maximum northward and 
southward excursions of the HCS (see |http://wso. Stanford. edu7Tilts. html |). At the moments shown 
in Fig. 4, the HCS tilt angle is 6.5° in the activity minimum period and 71.3° in the maximum 
period. 

A comparison between the tilt angle of the HCS resulting from our SFT model with the values 
derived by PFSS extrapolations of WSO maps of the surface magnetic field is shown in Fig. 5a. 
The results from our standard CSSS extrapolation (solid red line) do not significantly deviate from 
those obtained by a PFSS extrapolation of the same SFT results with R^s = 3.25Rq (dashed red 
line), suggesting that the calculated tilt angle does not sensitively depend on the extrapolation 
method (in contrast to the distribution of open flux). The other two curves in Fig. 5a represent 
PFSS extrapolations based on observed photospheric magnetic fields using two different boundary 
conditions for the photospheric magnetic field (see http://wso.stanford.edu/Tilts.html). Both the 
'classic' (green curve) and 'new' (blue curve) extrapolations provided by the WSO website assume 
that the field is potential between the photosphere and the source surface. The two models differ 
in the way the photospheric field observations are used as the inner boundary condition and in the 
height of the source surface. In the 'classic' case the surface observations are taken to correspond 
to the line-of-sight component of the potential field, and the source surface is assumed to be at 2.5 
Rq. In the 'new' model, the observed field is matched to the radial component of the potential 
field projected onto the line of sight: the assumption is that the field in the photosphere is purely 
radial and becomes potential only above the surface. Extrapolation for the 'new' case exist on the 
WSO website for source surface heights of 2.5 Rq and 3.25 Rq. We show the result for 3.25 Rq 
which is claimed to better match Ulysses data. 

The result from our model is consistent with both extrapolations. Figure 5b shows the dif- 
ference between the HCS tilt angle derived from the SFT model with CSSS extrapolation and the 
other three cases. The agreement with the WSO 'new' method is somewhat worse than that of 
the 'classic' method, especially around the activity minimum periods. For physical reasons, it is 
expected that the 'new' method assumi ng a purely radial photospheric field should be a better 



representation of the real solar situation (jWang &: Sheelevlll992l : iPetrie &: Patrikeeva II2009I ). Note 



however, that inferring the radial field strength from the measured line-of-sight component requires 
quite large corrections factors at high solar latitudes together with the larger uncertainties of the 
high-latitude data. This suggests that the tilt angles derived from observed synoptic maps during 
solar minimum periods should be considered with some caution. 
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3.4. The minimum of cycle 23 



The current solar minimum appears to be rather extended, with particularly low activity levels. 
There are in dications that the polar field strength is significantly lower t han during the two previous 
minima (e.g. lSvalgaard et al.ll2005l : ISchattenll2005l : ISchriiver k. Liull2008l ). A low polar field strength 
would be consistent with the small measured values of the near-Earth interplanetary magnetic field 
(and thus also the open fiux, see Fig. 3) and the relatively large tilt angle of the HCS inferred by 
the 'classic' PFSS extrapolation (blue curve in Fig. 5a). However, the uncertainties in the inferred 
values of the polar field are large and the HCS tilt angle determined with the (presumably more 
relevant) 'new' PFSS model using a radial-field photospheric boundary condition does not show 
unusually large minimum values. 

Our SFT model yields a polar field during the present minimum which is is not significantly 
weaker than that of the previous minimum (see Fig. 2). The open fiux determined from the model 
is consistent with the low 0MNI2 measurements during the current minimum; however, the model 
results tend to be too low during minima, so that the present agreement might well be fortuitous. 
The tilt angle of the HCS derived from the SFT model decreases to low values during the current 
minimum, which is consistent with the 'normal' polar field strength that the model yields. The 
tilt angles inferred from PFSS extrapolations based on observed photospheric field distributions 
provide a confusing picture (cf. Fig. 5a): while the 'classic' line-of-sight boundary condition leads 
rather large tilt angles, which would be in accordance with a weak polar field, the 'new' radial-field 
boundary condition predicts values not much different fror n those during previous cy cles. As the 
'new' method is considered to be physically more realistic (jPetrie &: Patrikeeva II2009I ). this would 
dilute the case for an unusually weak polar field during the present minimum. 

In the SFT model, the magnitude reached by the polar field in the second half of a cycle 
crucially depends on the tilt angle (with respect to the East- West direction) of the emerging bipolar 
magnetic regions during the cycle. We have assumed the same relationship, 7 = 0.15A, between tilt 
angle and emergence latitude for all cycles considered. However, the analysis of sunspot group data 
indicates that the factor of proportionality in this relation may actually vary from cycle to cycle 
(Dasi Espuig et al., in preparation). Systematically smaller tilt angles dur ing cycle 23 cou l d lead 
to a unusually weak polar field during the present minimum. Alternatively, ISchrijver &: Liul (120081 ) 
proposed an increased strength of the diverging meridional flow near the equator as an explanation 
for the weak polar field. 



4. Conclusions 

We have simulated the temporal evolution of the Sun's total open magnetic flux and the 
heliospheric current sheet (HCS) since 1976 by coupling a surface flux transport (SFT) model and 
the current sheet source surface (CSSS) extrapolation method, using the observed sunspot groups 
to provide the magnetic flux input for the model. We draw the following conclusion from our 
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results: 

1) The simulated open flux matches the 0MNI2 data quite well, except for systematically 
lower values in the ascending phase of the activity cycle. This is consistent with the view that the 
solar open flux is largely determined by the instantaneous photospheric sources. 

2) At the source surface, the magnetic fleld satisfles the 'split monopole' conflguration suggested 
by the Ulysses out-of-ecliptic measurements. 

3) The temporal variation of the tilt angle of the HCS from the SFT/CSSS model matches the 
values derived from potential field source sTirfacc(PFSS) extrapolations of the observed photospheric 
magnetic field. The best agreement is found for the 'classic' line-of-sight condition. 

4) The conditions during the present minimum period of cycle 23 as provided by the SFT/CSSS 
model are similar to those at previous minima: the polar field has about the same strength as that 
of cycle 22, the tilt angle of the HCS is small, and the open flux is roughly at the level of the last 
minimum. Since 2007, however, the HCS tilt angle has deviated significantly from the WSO PFSS 
values with line-of-sight boundary condition ('classic' case) and approached those with the 'new' 
radial boundary condition. These results may well be affected by a systematic variations of the 
sunspot group tilts with respect to the E-W direction from cycle to cycle, as indicated by recent 
analysis of sunspot observations. 

5) In spite of some deviations in detail, the overall agreement of the model results with ob- 
servationally inferred values of open flux and current sheet geometry is encouraging. It opens the 
possibility to extend the model backward in time by using the sunspot group record since 1874. 
This will be the topic of a subsequent paper. 

Acknowledgments: Y.-M. Wang kindly provided the observational datasets of the averaged 
unsigned photospheric field shown in Fig. 2. M. Lockwood kindly provided the kinematically 
corrected 0MNI2 data shown in Fig. 3. 
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Fig. 1. — Number of sunspot groups per month (upper panel) and time-latitude plot of the emerging 
sunspot groups (lower panel) that constitute the input sequence of BMRs to the SFT model . Data 
source: USAF/NOAA. 
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Fig. 2. — Upper panel: Averaged unsigned photospheric field as a function of time. The symbols 
show the observed data from the Wilcox and Mount Wilson observatories (averaged over Carrington 
rotations). The solid curve represents 27-day averages from the SFT simulation. Lower panel: The 
temporal evolution of south (solid curve) and north (dashed curve) polar fields (averages poleward 
of ±75° latitude). 
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Fig. 3. — Temporal evolution of the unsigned radial field at 1 AU. Shown are the results derived from 
the SFT/CSSS model (red curve) in comparison with near-Earth measurements (blue curve). The 
latter have been obtaine d by first averaging the (signed) 0MNI2 data over 1-day intervals to remove 
small-scale fluctuations (jLockwood et al.ll2006l ) and then carrying out a 3-month running average of 
the unsigned values. Th e green curve represeiit s yearly averages of the measurements after applying 
a kinematic correction (jLockwood et al.l 12003 ) to remove effects due to the longitudinal structure 
in the solar wind. Dashed and dotted vertical lines indicate the epochs of solar cycle maxima and 
minima, respectively. 
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Fig. 4. — Magnetic field distributions at solar minimum (left panels, in 1987.2) and solar maximum 
(right panels, in 2000.5) resulting from the coupled SFT/CSSS model. Maps of the radial field 
are shown at the solar surface (a,b), at the cusp surface at i?cs = l-Si?© (c,d), and at the source 
surface at Rss = WRq (e,f). Panels g and h show latitude profiles of the field strength at the source 
surface, taken at the longitude with the largest latitude excursion of the heliospheric current sheet 
(HCS). 
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Fig. 5. — (a) Temporal evolution of the HCS tilt angle. Solid red curve: SFT result with CSSS 
extrapolation (a = O.2i?0, i?cs=i.8R0,Rss=iORQ); dashed red curve: SFT result with PFSS extrapola- 
tion (i?ss = 3.25i?0); blue curve: PFSS extrapolation of WSO synoptic maps with line-of-sight field 
boundary condition ('classic' model, Rss = 2.5Rq); green curve: PFSS extrapolation of WSO syn- 
optic maps with radial field boundary condition ('new' model, i?ss = 3.25i?0). The data for the blue 
and green curves have been obtained from the WSO website (http:7/wso. stanford.edu/Tilts. html), 
(b) The difference between the HCS tilt angle from the SFT/CSSS and the three other methods 
shown in the same color scheme as used in (a). 



